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Abstract 

We discuss the hnal state interaction effects at high energies via a multi- 
channel N/D method. We find that the 2 by 2 charge-exchange final state 
interactions typically contribute an enhancement factor of a few times 10~^ 
in the B meson decay amplitudes, both for the real and the imaginary part. 

We also make some discussions on the elastic rescattering effects. 
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There have been increasing interests in recent years in studying the final state interaction 
(FSI) effects in B meson hadronic decays. In Ref. [0 , a Regge pole model was suggested 
to estimate the FSI effects in B decays and it was found that, for B —>■ DK decays the 
FSI effects are small in the sense that the corresponding partial wave S-matrix is close to 
unity. Especially the final state rescattering effects in charge-exchange processes are of order 
0(10^^). These results imply that for the color non-suppressed charged final states, a tree 
level calculation to the B decay amplitude can be a good approximation |^ . Using a similar 
method as in Ref. [|l|, the authors of Ref. |^ observed that, as a consequence of Pomeron 
exchanges the FSI effects do not vanish as the center of mass energy square s approaches 
infinity and therefore they suggest that there are sizable strong phases generated from FSIs. 
On the other hand, the authors of Ref. pointed out that in above estimates the real part 
contributions of the rescattering amplitudes were not considered. They use a dispersion 
relation to estimate also the real part contributions of the charge-exchange processes {B — >■ 
D-K^-KTi ) and find that Reggeon contributions (in charge-exchange rescattering processes) 
are one order of magnitude larger than those estimated in Ref. [00]. Since this result, if 
correct, leads to important phenomenological consequences [^, it is therefore worthwhile 
to re-examine the problem by carefully studying the FSI effects, including the real part, of 
charge-exchange rescattering processes, which is what this paper mainly devote to. We will 
also discuss the FSI effects of elastic rescatterings although there exist some uncertainties 
related to such rescatterings in the scheme presented in this paper. 

To correctly evaluate the FSI effects to the decay amplitudes, it is important to note that 
if the initial interactions are of short range the production amplitude and FSIs are actually 
factorized 0. Assuming (quasi) two-particle unitarity we will use the multi-channel N/D 
method to study the FSI effects. The inelasticity effects caused by those multi-particle 
states is presumably not large, since the pion multiplicity is rather low at s = m\ as the 
multiplicity increases only as log(s) (see Q and references therein) and at least part of 
the multi-particle final states can be classified as two-particle final states of resonances 
(accompanied by cascade resonance decays). 

2 



For a multi-channel system we can rewrite the partial wave T-matrix as, 

T = N/D , (1) 

where the bold face represents a matrix and N contains only the left-hand singularities and 
D the right-hand singularities. In the physical region, T and the production amplitude. A, 
satisfy the following unitarity relations, 

ImT = T*pT, ImA = A^T , (2) 

which hold when the center of mass energy square, s takes any value along the positive real 
axis above the lightest threshold in the complex s plane, if there is no anomalous threshold. 
We will come back to discuss the influence of the anomalous thresholds on high energy FSIs 
later in this paper and will argue that their effects are negligible. The unitarity relations 
lead to the following simultaneous integral equations for N and D, 

N=i/'''':'°"''..-.D^i-i/''":w'v. 

Tf Jl s' — s n JR s' — s 

where b is the discontinuity of T from the left-hand cut. Equation may be solved by 
iteration method. 

The FSI contribution to the production amplitude obeys a simple formula, 

A = Ao/D , (4) 

where A is the full decay amplitude and Aq is analytic in a region surrounding the positive 
real s axis. In the present physical situation, it implies Aq is a smooth function of s in a large 
domain surrounding s = m\ in the presence of inelasticity. It should be noticed that the 
dispersion integrals in Eq.(|^) for those components of N and D in which Pomeron exchanges 
involve should receive one subtraction. In potential scattering theory D is normalized to 
unity when s equals to infinity. Since the total cross sections of hadron-hadron scatterings 
are continuously rising as a function of s, it is impossible to have such a normalization 
condition for high energy hadron interactions. Aq in Eq. @) is the Born term in potential 
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scattering theory and be constant. In the present case there is an uncertainty in relating 
Ao to the bare amphtude. This is related to the multiplicative ambiguity in N/D method 
and the lacking of a natural normalization point in D. However, as will be seen below, this 
ambiguity can actually be avoided in calculating charge-exchange rescatterings. 

To extract the information on charge-exchange processes from Eq. we may rewrite 
D as, 

D = D - D% (5) 

where D contains only channels dominated by Pomeron exchanges and D'^ represents charge- 
exchange processes. Assuming D'^ to be a small quantity and that D does not deviate far 
from unity (these assumptions will be justified by the results given later), we can make an 
expansion in the denominator of Eq. in powers of D*^ and keep only the leading term to 
obtain, 

■^charge— exchange -^e^D /D , (6) 

where A^i is the physical amplitude in the limit of vanishing charge-exchange rescatterings: 

Ae, = Ao/D , (7) 

which can be a good approximation to the physical amplitudes for decay channels with 
large(larger) bare amplitude. Equation can be written in an explicit form after a little 
algebra, 

1 rp{s')T%s')V^{s')^^^, ^ 1 rp{s')T^{s')V%s')^^,^ 



^ch.—ex. 



n J s' — s ir J s' — s 



where s = m?^ and. 

It is straightforward to check that T)^ satisfies the following integral equation, 

^ ' 71 J {S" - S'){S" - S) ' ^ ^ 
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and Eqs. (||), (|D and (|^) define an integral equation for V^, provided T''', T'^ and are 
known. Tlie above equations are valid to the leading order in the expansion in powers of T*^ 
and to all orders of Pomeron exchanges. 

In order to calculate the FSI contributions in charge-exchange processes it is required 
to determine V^{s') through Eq. (|TD| ). It is a difficult task since for a process i ^ j under 
consideration there can be many channels (j' j) dominated by Pomeron exchange exist 
contributing to i j through i j' ^ j- To proceed let's look at the elastic process 
(j' = j) first and approximate the solution of Eq. (|l^) by first order iteration result, one 
obtains, 

^'^^ TT J {s" - S'){S" - S) ' ^ ^ 

where 

TP(s) = — (12) 

' 167rso(A + a'log(s/sO) -m«V2) ' ^ ' 

where (3^ is the Pomeron coupling to matter and A = 2.82 is the parameter characterizing 
the form-factor of the Pomeron coupling, also we set = 1 in above. From Eq. ([Tl|) we find 
that 'D^{s) is a smooth function of s and deviates very little from unity for any reasonable 
value of s. Even for s being as large as W boson mass square, one still has, 

p^;(M^) = 0.99 + 6.1 X lO^^i (13) 

in D-K case. The D matrix elements of inelastic diffractive channels are expected to be 
the same order as — 1 ~ 10~^ but with different phases. Large cancellation is natural 
when summing over all the intermediate states j' which ensures terms like T^^iT^^ in Eq. (^ 
remains to be second order comparing with T^-. Therefore we can approximate Eq. (|| ) by, 

A,n.-e.. = Ae,{- / P(flE^ds' + ^pT^} , (14) 

TT J S' — S 

which is, in simple words, obtained by neglecting all the diffractive scattering effects in 
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charge-exchange rescatterings-Q Equation (^4|) is the basic formula for numerical calculations 
in the following.^ 

In order to compare with previous work, we now focus on the process, B D^tt^ through 
a D^n^ intermediate state. The D^n^ D^n^ scattering amplitude is, 

Td+^^^dOttO = 2V2p , (15) 

where the p-Reggeon contribution is parameterized as 

P = m ( - ) . n ' (16) 

in the small negative t region. We take Sq — 2a'~^ and also neglect the t dependence of 
(3^, according to Ref. 0]. The Reggeon coupling adopted in Ref. Q is obtained by relating 
(5r to the on-shell vector meson coupling in the t-channel physical process which however 
over-estimates the magnitude comparing with the one directly extracted from high energy 
scattering data. We take the value used in Ref. 0] for numerical calculation and will come 
back to this difference later. 

The approximation made above will be self-consistent if the final results on FSI effects 
are small which is indeed the case as we find, 



6.9 + 4.4i) X 10-2 . (17) 



Similarly for the tt+tt- ti^ti^ rescattering, we have, 

^Tr+vr-^vrOTrO = 4p , (18) 

in which sq = <y'~^ is taken and the numerical result is. 



^ This equation was used in Ref. Q without any justification. 

^ For vrvr rescatterings, T)(.l{M^r) = 0.95 + 0.20i, which is not as close to unity as in the Dtt case. 
This is not surprising as the diffractive FSI effects are stronger for vrvr rescatterings. However, 
Eq. (|l^ should still be good enough to work with in an order of magnitude estimate. 



These results are significantly smaller in magnitude than the one obtained in Ref. [||. Still, 
a more careful fit to the high energy scattering data suggest that the p-Reggeon coupling is 
further suppressed by a factor about 2.5 P] comparing with the previously used. With this 
suppression factor the parameter is close to the value used in Ref. |l| 0. We therefore 
suggest that Reggeon contributions to FSIs are typically of order of a few times 10~^, both 
for the real parts and the imaginary parts. 

Of course Regge parameterization can not always be correct when decreasing s down 
to thresholds even though we know from old experiences that it may indeed work well in 
charge-exchange processes at rather low energies. To estimate the violation of the Regge 
parameterization which we call S, one has, 



where pT'^ refers to the averaged value of the magnitude of pT'^ in the range of s which 
characterizes the region where Regge parameterization fails and the last inequality comes 
from unitarity. Taking s = lOGeV"^ we see that 6 is at best a few percent, in a very 
conservative estimate. 

Strictly speaking, the present method only applies to a restricted class of T matrices 
which contain no anomalous threshold. An anomalous threshold comes from dynamical 
singularities which occurs only if a loosely bound composite system of hadronic constituents 



^ When extracting from high energy scattering data it is probably better not to include in the 
fit the nucleon-nucleon scattering data which contain strong SED violating contributions and can 
not be explained by a simple Regge pole model. Actually, the fact that meson-nucleon data looks 
better in favor of a simple Regge pole parameterization is the main reason to support the use of 
a simple Regge pole model for meson-meson scatterings. On the other hand, absorption effects, if 
exist, only reduce the low partial wave FSI effects. 
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is involved 0. The typical mass scale of an anomalous threshold is the scale of the masses 
of the light hadrons in the decay products [|1^]. From a simple dimensional analysis we are 
convinced that these effects, if appear, are further suppressed in powers of 1/s in high energy 
scatterings. 

In above we have discussed the effects of charge-exchange FSIs, we obtained the relation 
between the rescattering amplitude and the physical amplitude in the limit of vanishing 
charge-exchange rescatterings (the bare amplitude renormalized by diffractive rescatter- 
ings). The ambiguity caused by the normalization of D is absorbed into A^i- In discussing 
the elastic rescattering it is unavoidable to face such an ambiguity. As having been made 
clear in above discussions that the ambiguity comes from the fact that there is no natural 
normalization point for D, unlike the case in potential scattering theory where D is normal- 
ized to 1 at s = oo. The only thing can be done in the present formalism is to compare the 
FSI strength between two energy scales. Here we define the FSI strength at s (= m^) to be 
1/Da(s) where Da is normalized to unity at scale A, and A is chosen somewhere below W 
boson mass and much larger than m^.^ From the previous estimation on Vei{M^), recall 
that D{si) normalized at S2 is equal to the inverse of D{s2) normalized at Si, we conclude 
that the FSI effects remain small even for elastic rescatterings at s = m^, at least in the Dir 
case.0 Here we defined in some way (maybe rather artificial) the absolute value of the phase 
of the elastic channel, but we want to stress that what truly matters in a physical quantity 
(like the CP violation observables) is the difference between two elastic phases of different 
isospins which is of non-diffractive nature and should therefore be a small quantity. 



See the discussion given in Ref. 0. 

^ Of course, a more careful numerical study calls for a calculation beyond the first order iteration 
in determining function D, including inelastic effects which can only be done under various as- 
sumptions and approximations. The random S-matrix approach [|ll| can be helpful in extracting 
out the essential feature of a system with many inelastic channels. 
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To conclude we find that charge-exchange rescatterings in B decays contribute typically 
an enhancement factor of order of a few times 10"^ ~ where A is the Wolfenstein param- 
eter. For elastic rescatterings, the FSI effects most likely remain small. B decays contain 
rich phenomenology and an enhancement factor ~ A via charge-exchange rescatterings leads 
an important role |^. In this sense the result of this paper, if correct, is in disfavor of an 
experimental extraction to the rich physical phenomena. 
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